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Abstract 

The rate of energy dissipation in solutions of the body-forced 3-d incompressible Navier-Stokes 
equations is rigorously estimated with a focus on its dependence on the nature of the driving force. 
For square integrable body forces the high Reynolds number (low viscosity) upper bound on the 
dissipation is independent of the viscosity, consistent with the existence of a conventional turbulent 
energy cascade. On the other hand when the body force is not square integrable, i.e., when the 
Fourier spectrum of the force decays sufficiently slowly at high wavenumbers, there is significant 
direct driving at a broad range of spatial scales. Then the upper limit for the dissipation rate may 
diverge at high Reynolds numbers, consistent with recent experimental and computational studies 
of "fractal-forced" turbulence. 
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I. INTRODUCTION 



A fundamental principle of modern hydrodynamic turbulence theory is that nonlinear 
interactions between Fourier modes of the velocity field can transfer energy from directly- 
forced large spatial scales, through the so-called inertial range, down to a small dissipation 
length scale where viscosity effectively consumes kinetic energy and transforms it into heat. 
This turbulent cascade process has been intensively studied experimentally, numerically, and 
theoretically (at various levels of mathematical rigor) since the first half of the twentieth 
century. See, e.g., the book by Frisch ^| for an introduction and entry into the vast literature 
on this subject, which is still the focus of much current research. 

One profound consequence of the cascade mechanism is the so-called dissipative anomaly 
wherein a finite and non-vanishing residual energy dissipation persists in the singular limit 
of vanishing viscosity, i.e., in the infinite Reynolds number limit. This phenomenon is 
quantitatively described as Kolmogorov scaling of the energy dissipation, namely 

(3 = — = O(Re ) asRe^oo (1) 

where e is the total energy dissipation rate per unit mass, i is an integral (large) length scale 
in the flow characterizing the domain or a large scale in the forcing and flow, U is a turbulent 
velocity scale, and Re = Ui/v is the Reynolds number with v denoting the kinematic 
viscosity. Sreenivasan has collected together relevant data illustrating Kolmogorov scaling 



in experiments 



Moreover, given precise definitions 



and direct numerical simulations 
of all the quantities involved, this (3 ~ Re Kolmogorov scaling has been shown to be an 
upper bound for (weak) solutions of the 3-d incompressible Navier-Stokes equations driven 
by sufficiently smooth — in particular, square integrable — body forces 

While the cascade picture of turbulence requires that energy be predominantly injected 
in a relatively narrow range of spatial scales, some researchers have recently performed 
experimental and computational studies of fractal- generated turbulence. These are flows 
driven by spatially broadband fractal forces with certain scaling properites that inject energy 
directly at a wide range of scales — most notably at small scales that could otherwise only be 
excited by the cascade. Such forcing can impose a self-similar structure on the flow that is 
independent of the turbulent energy cascade. If such forcing can be achieved experimentally 
then one can observe, and in principle control, the balance between the energy that has been 
directly injected and the energy transfered by the nonlinear mode interactions. 



Indeed, Queiros-Conde and Vassilicos |8j performed experiments by forcing fluid past a 
fractal object, an obstacle that was structurally self-similar over several scales. Staicu et 
al P| experimentally measured energy spectra and structure functions in the wake of grids 
of fractal dimensions 2.05, 2.17, and 2.40 in a wind tunnel, concluding that "there may be 
a direct relation between the scaling properties of the fractal object and the turbulence it 
creates". This is more easily investigated in direct numerical simulations where details of 



the flow field is directly observable. 
Mazzi and Vassilicos 



1(3 ] performed direct numerical simulations of stationary homo- 
geneous and isotropic turbulence in a fluid in a 3-d periodic box of size I driven by a 
velocity-dependent fractal body force f(x, t) with Fourier components of the form: 

f k (*) = ^ vm*)i iki iu k (t)i; 11 (2) 

( , |k| > kp 

where Uk are the velocity field's Fourier components, fk = whenever Uk = 0, and uo = 0. 
The scaling exponent ( is intended to characterize the fractal properties of the stirrer or 
obstacle, and the maximum wavenumber kp is to be thought of as the inverse of the spatial 
size of the smallest parts of the fractal stirrer. Mazzi and Vassilicos used numerical values 
for which the fractal forcing extended down to scales ~ kp 1 on the order of the Kolmogorov 
dissipation length rj = (z/ 3 /e) 1 / 4 . They observed that the bulk energy dissipation rate did 
not exhibit Kolmogorov scaling (3 ~ Re , but rather (3 ~ Re 1 corresponding to e ~ XJ A jv. 



Biferale et al [ll|, [12j performed numerical simulations of the 3-d Navier-Stokes equations 
with a stochastic body force that was white-noise in time but with a power law spectrum of 
spatial scales ~ They investigated small scale turbulent fluctuations and concluded that 
the statistics displayed two distinct qualitative behaviors. When the spatial spectrum of 
the forcing decayed sufficiently fast, the small scale fluctuations were universal in the sense 
that they were independent of the details of the force spectrum. This regime corresponds 
to conventional cascade dynamics. When the spatial spectrum of the forcing decayed more 
slowly, however, the small scale fluctuations were "force-dominated" with the cascade being 
overwhelmed by the direct excitation from the driving. Interestingly, they reported that 
this transition occurs at a value £ = — | of the scaling exponent corresponding to the 
boundary between (spatially) square integrable and "rougher" forcing functions without 
square summable Fourier coefficients. 



In this paper we derive rigorous upper bounds on the bulk energy dissipation e in an 
incompressible Newtonian fluid driven by a variety of body forces including forces that are 
not square integrable. This work generalizes the previous analysis for square integrable 
body forces 4 , Q, E | to include fractal forces that drive the flow directly at a broad range of 
scales. In accord with the findings of Biferale et al we find that the case of square integrable 
forcing is a borderline situation: (3 < Re when the body forces are square integrable (or 
smoother), but the estimates increase for rougher driving so that the dissipation coefficient (3 
may increase asymptotically as Re — > oo. For the roughest forcing functions that make sense 
mathematically, i.e., forcing functions with Fourier coefficients satisfying ^ k ^~ 2 |ik| 2 < <x>, 
we find that (3 < Re 1 , the scaling observed by Mazzi et al . 

The rest of this paper is organized as follows. The following Section II lays out the 
mathematical setting for the analysis and gives the definitions of the physically relevant 
quantities of interest. In Section III we study the case of time-independent body forces, 
and the subsequent Section IV deals with velocity-dependent forces like (j2j). The concluding 
Section V contains a brief summary and some closing remarks. For completeness and to 
make the paper self-contained, we include some mathematical details in an appendix. 



II. STATEMENT OF THE PROBLEM AND DEFINITIONS 

Consider the incompressible 3-d Navier-Stokes equations on a periodic domain x £ [0, £] 3 : 

d t u + (u- V)u + Vp = z/Au + f (3) 

where u(x, t) is the divergence-free velocity field, p(x, t) is the pressure, f (x, t) is the applied 
body-force, v > is the kinematic viscosity, and u| t=0 = u (x) is the initial condition. We 
will take the body force to be a specified time independent (divergence-free) function f(x) 
in Section III, or given by a velocity-dependent expression like in Section IV. We write 
Fourier expansions as 

u(x,t) = Vu k (t)e ix - k where u k (t) = ^ I e" ix - k u(x, t)d 3 x. (4) 

Without loss of generality, in every case we will take the applied body force and initial data 
to have spatial mean zero so that k = |k| > in all sums. 
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A field u(x, t) g H a if ||u(-, t)||#a < oo where we define the Sobolev norms || • by 

IM,*)ll2r- = ^(^) 2 "|u k (t)| 2 = I [\(-fAyu(x,t)\Wx. (5) 

The index a can be positive or negative and the function spaces H a are nested according to 
H a C H a for a > a'. The case a = corresponds to the usual L 2 norm (with the volume 
normalization) and we write 

||u||#o = ||u|| L 2 = ||u|| (6) 

"Fractal" forces are defined as those with power- law Fourier coefficients, |fk| = Cfr. For 
such a function to belong to the Sobelov space H a its exponent must satistfy £ < —a — |, 
i.e., the Fourier coefficients must decay as |fk| < k~ a ~^~ 8 for some 5 > 0. 

We define time averages of functions g(t) according to 

g= lim / g{t)dt (7) 
T ^°° Jo 

and for simplicity in this paper we presume that this limit exists for all quantities of interest. 
The bulk (volume and time) average of a function h(x, t) is denoted by 



(h) = i fc(x, -)^- (8) 
The root means square velocity U of a solution u(x, t) of the Navier-Stokes equations is 

f/=<|u| 2 ) 1/2 =hF 1/2 , (9) 
and the bulk energy dissipation rate (per unit mass) is defined by 



e = <HVu| 2 > = ^||u||^. (10) 

When a solution u(x, t) satisfies the energy equality (i.e., when the energy is absolutely 
continuous, which holds for every regular solution), the energy dissipation rate satisfies 

e=(f-u). (11) 

That is, the power supplied by the driving force is balanced by the viscous dissipation. 

Weak solutions to these 3-d Navier-Stokes equations exist for f g H~ l , and then in 
general the relation in ()11|) is only an inequality, i.e., e < (f • u) ^ 
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1511 . This fact does 



not affect our results, however, because we will just derive upper limits on e. Moreover, 



the assumption of the existence of the long times averages is not necessary if the limit is 
replaced by limsupj^^. With that definition the estimates we derive are fully applicable 
to weak solutions. 

Using the definitions above, the Reynolds number is identified Re — Ui/v and the dis- 
sipation coefficient as (3 — e£/U 3 . In the scenario described here both Re (or U) and j3 
(or e) are formally "emergent" quantities, not directly controllable but determined rather 
as functions of i and v and functionals of uo and f . These bulk averaged quantities gener- 
ally depend on u , but the relationships derived below are uniformly valid for all solutions 
regardless of initial data so we will drop any further reference to them. In practice one 
assumes that the parameters of the force, e.g., its amplitude, can be tuned to achieve any 
desired Reynolds number. Then f3 may be considered a function of Re. The overdamped 
highly viscous limit is Re — > and the vanishing viscosity limit is explored as Re —>■ oo. 

Some very general statements can be made for the overdamped limit. Poincare's inequal- 
ity implies that 

Aix 2 i> , „ 

e > (12) 

so for any forcing 

4tt 2 

P>^. (13) 

This Reynolds number scaling is sharp: as will be seen below, for a wide variety of forces 
there exists a constant c > 47r 2 (generally depending on the details of the forcing) such that 

< 4~ as Re -> 0. (14) 
Re 

This scaling, (3 ~ Re -1 , is characteristic of large scale laminar flows where the typical rate 
of strain is proportional to U jl and the typical stress is proportional to vU /£. 

For higher Reynolds numbers the lower estimate in (|13|) can generally not be improved. 
That is, at arbitrarily high Re there are forces that can sustain the flow with (3 ~ Re -1 . 
Those flows — which may be unstable — are necessarily sufficiently laminar to exclude any 
characterization as being turbulent. The upper bound on (3, however, necessarily increases 
above Re -1 as Re — > 00. For turbulent flows with an effective energy cascade the dissipation 
becomes independent of v as Re — > 00, i.e., (3 ~ Re , as evidenced by experiments and direct 
numerical simulations. But for sufficiently broadband forcing (3 may increase indefinitely in 
this limit, and the task of the next two sections is to place rigorous upper bounds on (3 as 
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a function of Re for flows driven by fractal forces. 



III. STEADY H~ a BODY FORCES 



In this section we generalize the approach introduced by Doering &: Foias 



that was inspired by previous work of Foias and coworkers 



16 



13 Q 



a 



-an approach 



-to cases where the 

time independent force f(x) G H~ a with a G [0, 1]. For a < the force f G L 2 and the 
(3 < Re upper bound, corresponding to the usual energy cascade, is effective |4|. We do 
not consider values of a > 1, for then even weak solutions of the Navier-Stokes equations 
are not known to exist. While the analysis in this section is not restricted to strictly fractal 
forces, the results apply nevertheless to those with power-law Fourier coefficients |fk| ~ 
where ( = a — | — S for any 5 > 0. 
Write the steady body force as 



f (x) = F^rSc), 



(15) 



where F is the amplitude of the force, the H a norm of f, and the "shape" function is a 
dimensionless divergence-free field on the unit 3-torus normalized according to 



U\\ H - a = i. 



(16) 



Using Cauchy-Schwarz and the interpolation inequality ()A6j) with s = a, r = 1, t 
estimate 



0, we 



1 



f • u dx 



£ 3 



< $^|fk| Kl < 

k 



1/2 



1/2 



|f||jr—||u||tf«. = F\\u\\ Ha < FlIull^iHull 1 - . 



Then taking time average and applying Holder's inequality, 



Note that \_ E [0, |] so Jensen's inequality (|A2|) ensures that the last term in the last 



2-a 

line above is bounded by U l ~ a . Hence 

2a a 2 2(l-g) 

e < £ 2 - Q z/ Z7 2-a . (17) 

On the other hand we can also estimate F from above independently in terms of U, v 
and t. Multiply the Navier-Stokes equation Q by a sufficiently smooth time- independent, 
divergence-free function , 0(£~ 1 x) on the unit 3-torus satisfying ((f) • if)) > 0. (It's easy to 
produce such fields if), for example as a finite Fourier mode Galerkin truncation of (p.) 
Integrating by parts, taking time averages, and applying Holder and Cauchy-Schwarz, 

rV2 



F((j>-i/)) = -(u-(Vif))-u)-D(u-Aif))< ||V^||l-||u|| 2 + z/||AV||||u|| 2 ' . (18) 
Hence 



F< 1 



((f). if)) £ 1 "-"--£2 j > ( 19 ) 

where V = £V is the dimensionless gradient on the unit 3-torus . Plugging this estimate for 
F into the bound (JTTj) for e we deduce 

P < Re^ (Ci + C 2 Re" 1 ) ^ , (20) 

where the coefficients Cj depend only on the shape function (f) and the multiplier function 
if) — but not on the parameters of the problem, i.e., the force strength F, the viscosity v, or 
the outer length scale L Specifically, 

C| J'^»^ and Q= MLjWk. (21) 
For i?e> 1 the upper bound (ffiij) scales 

/? < i?e^ (22) 

where the exponent £ [0, 1]. If a = 0, i.e., when the force f G L 2 , we recover the classical 
estimate corresponding to Kolmogorov scaling 

P < 1 (23) 

r-i 

that holds as well when o<0 j|. In the other borderline case a = 1, 

P < Re- (24) 
8 



And as advertised, when Re <C 1 the overdamped laminar scaling 

P < Re' 1 (25) 

emerges for all a < 1. 

IV. A TIME DEPENDENT FRACTAL FORCE 







Following Mazzi & Vassilicos |H| , consider a fractal forcing function of the form: 

f k (t) = F(£|k|)^f^ + i^i) , (26) 

VjUkl |k||u k |y 

where F is the strength coefficient, ( G [— |, — |] and 5 > and i*k = whenever Uk = 0. The 
Navier-Stokes equations (j3J) driven by this velocity-dependent time-varying force constitute 
an autonomous system. We assume initial data u (x) ^ 0, that a (statistically) steady flow 
is subsequently sustained, and that for t > each |uk(£)| = only on a measure- zero set of 
times. The scaling exponent £ = — | corresponds to the case where the forcing is in H^ 1 at 
each instant of time for all 5 > 0, while £ = — | (or less) is L 2 (or smoother) forcing. 
Start by writing 



e = (f • u) = f k • ui 



F £p^(« c+3/2 i^i- ( 27 ) 



Applying the Cauchy-Schwarz inequality, 

-1/2 



e^FC^^k^C+aiukl^ . (28) 



where 

1/2 



c 



k (W) J 

Note that the (3-d) sum defining C converges iff 5 > 0. Indeed, C = C(5 -1 / 2 ) as 5 — > 0. 
Holder's inequality then implies 

-C/2-1/4 / x C/2+3/4 



u k | 2 (30) 



t p, \ C/2+3/4 



CFU-^ 1 ' 2 ( ^- j . (31) 
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Solving for e, 

2C+3 

e<C—F—U—^-y Q . (32) 

Now the challenge is to eliminate F in favor of U, i and v. 

To derive an upper bound on F we will estimate the bulk average of the Q dotted into 
the time-dependent test function i/>(x, t) with the Fourier coefficients: 

Mt) = 7^|k|)- 4 - 5 ' (33) 

l u k| 

for |u k | 7^ 0, with ip-^ = when |u k | = 0, and 5' > 0. We consider the resulting terms one 
by one. 

First note that the pressure term (if) ■ Vp) = since V • ip = 0. The advection term is 
estimated 

| (V ■ (u ■ Vu)> | = | (u ■ (VV>) • u) | < ||V^|U-c <|u| 2 > (34) 

where 

llv^iUo. < r^iki |y, k (t)| < r l J2( £ \ k \)~ 3 ~ s ' = t ( 35 ) 
k k 1 

and the pure number C is finite for all 5' > 0; C = 0(5'^ 1 ) as 5' — > 0. The force term is 

(f • V) = F ^(£|k|)^ 4 - <5 - <5 ' = C"F (36) 

k 

where the sum for C" converges uniformly for all non-negative 5 and 5' and all ( < —1/2. 
By asserting equality in (J3H|) above we use the assumption that |u k | > for almost all t > 0. 
Next, the viscous term is 

■ ^u) = v (AiP • u) < [/ = C" ^ (37) 

where ^ 

C "' = fe £ l k l _4_25 'j ( 38 ) 
is uniformly bounded for all 5' > 0. Finally, observe that 



(^•9 t u) = ^|u k |(£|k|)-^'. (39) 

k 

The time average of the time derivative of a quantity vanishes if the quantity is uniformly 
bounded in time. Because 

$> k (t)i(£ikir 4 -*< (5> k (oi 2> ) ' (j2mr 8 - 2S ') ' m 
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where |uk(t)| 2 = ||u(-,t)|| 2 is uniformly bounded in time for these H^ 1 (or smoother) 
forces, the sum above converges for all 5' > and we conclude that ■ <9 t u) = 0. 
Hence the bulk average of if) dotted into the Navier- Stokes equations yields 

F< cv ™ 

~ C" £ C" P 1 ; 

with absolute constants C" and C", and C depending only on 5' > 0. Inserting into (|52*|). 

3+2C , -i \ 4 

(3 < Re~ (ci + c 2 Re~ l ) ^ (42) 

where 

c i = "^T and c 2 = "^T- ( 43 ) 
As before, when Re — > 0, this result produces the laminar scaling 

P < (44) 

for all relevant values of the force's scaling exponent. When Re — > oo, however, the dissipa- 
tion may be as large as 

P < Re^ (45) 

with exponent < |^|| < 1 as ( varies from — | to — ~. It is worthwhile noting that the 
coefficients c\ and c 2 depend on 5 > (and 5' > 0, introduced for convenience) — but not at 
all on the force parameters F and £ or on v or I — and that the coefficients c\(S) and 02(6) 
diverge as S — > because C(S) defined in diverges as 6 — > 0. 

V. SUMMARY & DISCUSSION 

In this paper we generalized the analysis that was previously employed for square inte- 
grable (or smoother) steady forces (^J and velocity-dependent forces 7| to derive bounds on 
the energy dissipation in the case of broad-band and fr act ally-forced flow described by the 
incompressible 3-d Navier-Stokes equations. When a steady body-force f(x) e H~ a with 
a G [0, 1], we showed that the dimensionless dissipation factor /3(Re) is limited according to 

An 2 Re' 1 < [3 < Re^ (d + C 2 Re a ~ 2 ) ^ . (46) 
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For velocity-dependent fractal forces of the form (|2l)]) with |fk| ~ k^ s , ( G [— §, — 5], and 
5 > 0, we deduced that 

3+2C 4 

47r 2 i?e- 1 < (3 < Re— ( Cl {5) + c 2 (5)Re- 1 ) ~ . (47) 

These scalings are sharp as Re — *> 0, displaying the laminar behavior (3 ~ Re" 1 . As 
Re — > 00, both upper estimates are (3 ~ Re for square integrable forcing, i.e., a = and 
£ = — |. And in the extreme limits a = 1 and £ = — |, both estimates give (3 ~ ite 1 . We 
remark that the scalings in (j4*U|) and (|4T|) are clearly consistent with each other when it is 
recognized that forces with |fk| ~ k^~ s are in H~ a when £ = a — |. 

In terms of dimensional physical quantities, we have estimated the energy dissipation rate 
(per unit mass) e in terms of the rms velocity U, I and v. Laminar dissipation corresponds 
to e ~ uU 2 /£ 2 while the turbulent cascade is characterized by e ~ U 3 /£ and the roughest 
fractal forces may allows e ~ U 4 /u. But for a specified form of the bodv force it is natural 
to consider e and U as functions of the forcing amplitude F, £ and v [l9(. When the force 
is specified, rather than the rms velocity, it is well known (and easy to show) that the 
Stokes flow driven by the given force sets an upper limit for the dissipation rate; any other 
flow necessarily dissipates less energy. In terms of the F, i and u, the maximal Stokes flow 
dissipation is e ~ F 2 £ 2 / v which may be interpreted as a laminar flow bound. It is interesting 
to note that in the extreme limits of H~ l forcing in (|T7|) and £ = — | in ()32|) . the scaling in 
this laminar upper limit is reproduced explicitly. 
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APPENDIX A: INEQUALITIES 



For convenience, in this appendix we collect the mathematical estimates used here: 

(a) Jensen's inequality: If the real- valued function of a real variable 9(x) is convex, then 
for each real- valued function g 



e({g))<(8og) 



(Al) 



where (•) stands for averaging. In particular, for any nonnegative function g and any 
real number p E [0, 1], 

(g p ) < (g) p • (A2) 



(b) Holder's inequality: 

</>(x) VK X ) dx 



< ( / |</>(x)('"dx] Y/"|^(x)|Mx N 



(A3) 



valid for all 4> E L p and tp E L q , where p and q > 1 and ^ + - = 1. For an V sequence 
(a&) and an / 9 sequence (where p and g are related as above) the discrete analogue 
of (lA"3j) reads 



Vp 



1/9 



< 



An important case of (jA4|) (for p = q = 2) is the Cauchy-Schwarz inequality 



(A4) 



(A5) 



(c) Interpolation inequalities between Sobolev spaces: Let < r < s < t and u E H l . 
Note the algebraic identities 

t — s s — r 



t — s s — r 
r h t = s, 



+ 



I. 



t — r t — r t — r t — r 

These interpolation estimates are the result of applying Holder's inequality (jA4|) in 
Fourier space: 

hi* = ^w 2 iu k r = ^[(^) 2 nu k r]^[(^) 2t |u k | 2 ]^ 



< 



2 t-s 2 s^r 

lullr'^Hulli t_r . 



S — 7' 

t—r 



(A6) 
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